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Paramagnetic impurities in a quantum spin-liquid can result in Kondo effects with highly unusual 
properties. We have studied the effect of locally exchange-coupling a paramagnetic impurity with the 
spin-^ honeycomb Kitaev model in its gapless spin-liquid phase. The (impurity) scaling equations 
are found to be insensitive to the sign of the coupling. The weak and strong coupling fixed points are 
stable, with the latter corresponding to a noninteracting vacancy and an interacting, spin-1 defect for 
the antiferromagnetic and ferromagnetic cases respectively. The ground state in the strong coupling 
limit in both cases has a nontrivial topology associated with a finite Z 2 flux at the impurity site. For 
the antiferromagnetic case, this result can be obtained straightforwardly owing to the integrability 
of the Kitaev model with a vacancy. The strong-coupling limit of the ferromagnetic case is however 
nonintegrable, and we address this problem through exact-diagonalization calculations with finite 
Kitaev fragments. Our exact diagonalization calculations indicate that that the weak to strong 
coupling transition and the topological phase transition occur rather close to each other and are 
possibly coincident. We also find an intriguing similarity between the magnetic response of the 
defect and the impurity susceptibility in the two-channel Kondo problem. 


I. INTRODUCTION 

A study of disorder in condensed matter systems is use¬ 
ful from two perspectives. Disorder is inherent in most 
condensed matter systems and often has profound effects 
on their properties. Incorporation of small amounts of 
paramagnetic impurities in a metallic host can result in 
the Kondo effect which gives the well-known logarithmic 
temperature dependence of the resistivity upon cooling, 
and eventually crosses over to a Fermi-liquid regime with 
a characteristic low energy scale, the Kondo tempera¬ 
ture. Conversely, impurities at low concentrations can 
act as a probe providing specific signatures of the envi¬ 
ronment they exist in. From the latter perspective, the 
Kondo effect is a set of signatures of certain low-energy 
excitations of the host lacking long-range magnetic or¬ 
der. For instance, exotic Kondo effects are known to 
arise in itinerant electron magnets near criticality^^— and 
in insulating quantum spin-liquid systemsA.^ owing to 
the paramagnons and spinonic excitations respectively. 

A study of impurity effects in the spin-^ honeycomb 
Kitaev model— is very appealing in this context. This 
Kitaev model is integrable and the ground state can 
be either a gapless or gapped quantum (Z 2 ) spin-liquid 
with extremely short-ranged spin correlations— The ele¬ 
mentary excitations are not spin-1 bosons that one typ¬ 
ically expects for magnetic systems in two dimensions 
and higher, but emergent dispersing Majorana fermions 
(spinons) and Z 2 vortices (tt— flux excitations associated 
with spins at the vertices of the hexagonal plaquettes) 
which in the gapless phase are known to be non-Abelian 
anyons.— Experimental realization looks increasingly im¬ 
minent with several interesting proposals to realize Ki¬ 
taev physics in two-dimensional quantum-compass ma¬ 


terials such as the alkali iridates^ and ruthenium trichlo¬ 
ride;^ and independently in cold-atom optical lattices-^ 

Introducing a paramagnetic impurity into the model 
through local exchange-coupling of the impnrity spin 
with a host (Kitaev) spin results in a highly unusual 
Kondo effectii owing to the peculiar elementary exci¬ 
tations in the host. For an 5" = 1/2 Kitaev model with 
an energy scale J coupled locally to a spin-^ impurity, 
the perturbative scaling equations for the impurity cou¬ 
pling K turn out to be independent of its sign, with an 
intermediate coupling unstable fixed point |Ar| ~ J/S 
separating weak and strong coupling regimes. Such scal¬ 
ing differs qualitatively from the Kondo effect in met¬ 
als (and graphene/^) where a nontrivial effect is seen 
only for antiferromagnetic coupling, but is similar to the 
Kondo scaling reported for paramagnetic impurities in 
certain pseudogapped bosonic spin-liquids— The distin¬ 
guishing feature of the Kitaev-Kondo problem is that the 
weak and strong coupling limits correspond to different 
topologies of the ground stateJ^ 

Despite the insensitivity of the scaling equations to the 
sign of impurity coupling, the strong coupling limits for 
K > 0 and K < 0 are very different physically. In the 
antiferromagnetic case (K > 0), the strong-coupling limit 
for an S' = 1/2 impurity spin corresponds to a spin singlet 
at the impurity site - equivalent to the Kitaev model 
with a missing site, which is an integrable model. In the 
ferromagnetic case, the strong-coupling limit corresponds 
to a non-integrable problem where one of the sites has 
S = 1, while the rest have S = 1/2. 

The problem of missing sites (spinless vacancies) in 
the Kitaev model has received much attention in recent 
times. It was independently reported in Ref. [ll| and 
Ref. [ 1 ^ that the ground state of the Kitaev model with 
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a missing site is associated with a finite Z 2 flux through 
the defect. That would not be the case, for example in 
graphene, where although the Dirac fermions have the 
same dispersion as the emergent Majorana fermions in 
the Kitaev model, the phases of the intersite hopping 
matrix elements in graphene are identical for every bond 
and do not change upon the creation of defects. In con¬ 
trast, the phases of the intersite hopping elements of the 
Majorana fermions in the Kitaev model are a degree of 
freedom and can take values 0 or tt. For the missing site 
problem, the magnetic susceptibility is predicte d^^d^ to 
have logarithmic singularities both as a function of the 
magnetic field as well as the temperature. Some of the 
singularities in magnetic susceptibility reported in Refs. 

an dM are reminiscent of the two-channel Kondo prob¬ 
lem, and we shall later discuss a connection between 
such singularities and the presence of bound, zero en¬ 
ergy Majorana fermions in the Kitaev model with a va- 
cancyii as well as in the two-channel Kondo modelJ^ 
Vacancy induced spin textures have also been studied by 
exact diagonalizatioiii^ of finite clusters of up to 24 spins 
described by more general (and nonintegrable) Kitaev- 
Heisenberg models in the presence of a small magnetic 
field. In Ref. [13, it was reported that a vacancy induces 
longer ranged spin-spin correlations that extend beyond 
the single bond correlations one has in the defect-free 
Kitaev modeh§ The low-energy properties of the Kitaev 
model with a random and dilute concentration of va¬ 
cancies are also quite interesting. Such rare but locally 
singular perturbations are predicted to result in qualita¬ 
tively different low-energy properties compared to that 
expected for Gaussian white noise type disorder^ 

In contrast to the understanding we currently have on 
the effects of single and multiple vacancies in the Kitaev 
model, much less is known about the effect of spinful de¬ 
fects where the defect site has a different spin from the 
host sites. Part of the reason is that while the vacancy 
problem is integrable and affords a simplification of a dif¬ 
ficult problem where interactions and disorder are both 
present, the problem with an S' = I defect cannot be re¬ 
duced to a noninteracting one. Some progress was made 
in Ref. [ll| where it was explicitly demostrated that the 
ground states of both the vacancy as well as S = I defect 
problems have a two-fold degeneracy. However, it could 
not be established whether the S = I defect was also 
associated with a finite Z 2 flux that is the case when a 
vacancy is present. It was also not demonstrated whether 
the strong coupling fixed points were indeed stable. The 
stability is an important issue, for otherwise one would 
expect new intermediate coupling stable fixed points and 
not only would our understanding of the Kondo effect in 
the Kitaev model be incomplete, but also the paramag¬ 
netic impurity route for generating vacancies and spinful 
defects would no longer be appropriate. The latter is¬ 
sue is of interest from a practical point of view too since 
it would make it possible to generate nonabelian anyons 
conveniently using a spin-polarized STM tip to bind a Ki¬ 
taev spin ferromagnetically or antiferromagnetically. It 


is also not clear whether the topological transition and 
the magnetic transitions are coincident or occur at the 
same value of the impurity coupling strength. 


In this paper we address these open questions and 
make the following new findings. We demonstrate the 
stability of the strong coupling limit demonstrated for 
the ferromagnetic and antiferromagnetic cases which im¬ 
plies stability of the spin-0 vacancies and spin-1 defects 
created through this route. We perform exact diagonal- 
ization calculations for a finite fragment of the Kitaev 
model coupled to a paramagnetic impurity and show that 
while for weak impurity coupling, the ground state cor¬ 
responds to zero Z 2 flux at the impurity site, for strong 
coupling, the ground state has a finite flux irrespective of 
the sign of impurity coupling. For the value of impurity 
coupling that corresponds to this topological transition, 
we also observe the total spin at the impurity site going 
to zero or one depending on the sign of coupling - this 
establishes that the the topological as well as the weak 
coupling to strong coupling transitions occur very close 
to each other and possibly at the same point. As a corol¬ 
lary, we find that the ground state with a spin-1 defect 
corresponds to a finite flux at the defect site. Finally, we 
report an intriguing connection between the susceptibili¬ 
ties of a vacancy in the Kitaev model and of the magnetic 
impurity in a two-channel Kondo model. 


The rest of the paper is organized as follows. SecHIl 
provides a brief introduction to the honeycomb spin-^ 
Kitaev model. In Sec lIIIl we study the effect of coupling 
an external paramagnetic impurity to the Kitaev host 
through a local exchange coupling. A poor man’s scal¬ 
ing analysis of the impurity coupling reveals an unsta¬ 
ble fixed point separating the weak and strong-coupling 
regimes. The stability if the strong coupling fixed point is 
demonstrated for both ferromagnetic and antiferromag¬ 
netic impurity couplings. It is also shown that the strong 
and weak coupling limits correspond to different topol¬ 
ogy of the ground state. New conserved quantities are 
identified which are composite operators of an impurity 
spin component and two Kitaev flux operators. Sec lIVI 
contains the result of exact diagonalization calculations 
of finite Kitaev fragments coupled to external spins. The 
key findings in this section are (a)establishing that a fi¬ 
nite Z 2 flux is associated with the ground state of the 
spin—i Kitaev model with a spin—I defect - just as in 
the case of a vacancy, and (b) the magnetic (Kondo) 
and topological transitions occur very close to each other 
and are possibly coincident. In Sec|V]we discuss the in¬ 
triguing parallels between the low-temperature magnetic 
response of the Kitaev model with a missing site and the 
two-channel Kondo model. SecjVl] contains a discussion 
of the results and possible future directions. 
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Figure 1. (a) Schematic of a fragment of the Kitaev lattice 
showing the A and B sites and the x, y and z types of bonds, 
(b) Figure showing the reciprocal lattice vectors for the A sub¬ 
lattice. The Dirac point for the massless Majorana fermions 
is denoted hy kp and momentum summations are over the 
(shaded) half Brillouin zone. 


II. THE SPIN-1/2 KITAEV MODEL ON THE 
HONEYCOMB LATTICE 

Kitaev’s spin-i honeycomb lattice model for a quan¬ 
tum spin liquid is a model of direction dependent nearest 
neighbour exhchage interactions on a honeycomb lattice.— 
The Hamiltonian for this model is given by 

Ho = -J^Y1 (1) 

x-links y-links 2-links 

where the three bonds at each site (see FigH]) are la¬ 
beled as X, y and z. The model is exactly solavable— 
As was shown by Kitaev, the flux operators Wp = 
defined for each elementary plaquette p 
are conserved (FiglT]), with eigenvalues ±1, and form a 
set of commuting observables. The Kitaev spins can be 
represented in terms of Majorana fermions , Ci 

as erf = ibfci- This representation spans a larger Fock 
space, and we restrict to the physical Hilbert space of 
the spins by choosing the gauged Di = ibfbfb^a = 1. 

For each a—type bond, = ibfb'j is also conserved 
and the flux operators can be written as a product of 
Uij’s on the plaquette H^jj^gpiaq ity. The ground state 
manifold corresponds to a vortex-free state where all Wi 
are equal. In the vortex-free state, we can fix all Uij = I 
(corresponds to Wp = I) and the Hamiltonian can be 
written as a tight-binding model of noninteracting Majo¬ 
rana fermions. The reduced Hamiltonian for this ground 
state manifold is given by Hg = i ^JkCjCk, where 
Ajk = if k are neighboring sites on an a—bond 

and zero otherwise. The excited states (with finite vor- 
ticity) are separated from the ground state manifolds by 
a gap of order Jq,. 

The free Majorana fermion hopping Hamiltonian can 
be diagonalized in momentum space by defining the Bra- 
vais lattice with a two-point basis (FiglT]). In momentum 
space, 

Ho £o!(*l)®q,a®q.«J (2) 

q>0,Q: 


withe«(q) = ±|/(q)|, /(q) = 2f(+ 
Jz), where ag is the nearest neighbor spin distance. The 
eigenstates are a^^g = Cq_A + and Oq^ = 

Cq,A — with a{q) being the phase of /(g). 

The Kitaev model has gapless excitations for a region 
of parameter space where J^s satisfy the triangle in¬ 
equalities \Jx\ + \Jy \ > \Jz\ etc. and a gapped spectrum 
outside this parameter regime. The gapless phase has a 
point Fermi surface where e(kir) = 0 and e(q) has a lin¬ 
ear dispersion around kp (Fig[T]). For simplicity, we will 
assume Jx = Jy = Jz = J for further analysis. For this 
case, the Fermi points are at (±47r/3-\/3ao, 0). 

The ground state of the Kitaev model is a quantum 
spin liquid with only nearest neighbor spin-spin corre¬ 
lations.— On an a—bond, only (crfcr“) is non zero and 
other two spin correlations are zero. Four spin bond- 
bond correlations are long-ranged with power-law decay 
in the gapless phase of the Kitaev model. 


HI. TOPOLOGICAL KONDO EFFECT 


Consider a spin S impurity locally exchange-coupled 
to a host (Kitaev) spin at an A site (r = 0): 


Vk =* ^ E ^5“6“(dq,A + <,a) 

a qGHBZ,a 

E Q“5“6“(aq./3 + <^). (3) 

q^HBZ,a,p 


We perform a poor man’s scaling analysi o^^d^ for 
the Kondo coupling K to study the screening of the 
impurity spin by the host excitations. To study the 
system properties at low temperatures, we can com¬ 
pute the effective Hamiltonian for a reduced bandwidth 
for the fermionic excitations {—D -\- 6D to D — 6D) 
by integrating out the excitations in the band edges 
{{—D to — D + SD) and {D to D — SD)). This process 
is successively repeated to get a scaling law for the cou¬ 
pling constants in the Hamiltonian. We consider the 
Lippmann-Schwinger expansion for the T—matrix ele¬ 
ment, {il,b^\K^Sf^b^Ca. A|T2 + (q, Q;)). making a perturba¬ 
tion expansion T = t( 2) _(-... in increasing powers 

of K and following its variation as a function of the de¬ 
crease of the bandwidth {—D,D), we find that the first 
correction to the bare T—matrix comes from two third 
order terms (see Fig. HI). The contribution from on-site 
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scattering (Figl2^) is 

= {ny\VKG+{E)VKG+{E)VK\^ + {ci,a)) 




N 

iD-6D)<\ 



X G+(e) b^ia 


{Qh^iS^r{b^\bl{a^. 

q'.5+a|,,_5)G+(e) b^c 


+ aq//,5') 
q,Q: I (Qi O^)) 




N 


Y,{Q^f{SPf ^flqMaqM E-{Ho- e,M) 
q'J 


- £6 


~ -2Q^S^ 


p{D)a^\SD\ 

E-D 


— Y^iQP)\sPf 


E-J 


(4) 




Figure 2. Diagrams contributing to the scaling of Kondo cou¬ 
pling K°^. 


Here p{D) is the density of states at the band edge, a is 
the lattice constant and Go{E) = {E — Hq + iS)~^. 

Similarly, the contribution from Fig. [2jb) is 
y(3),f, ^ 0^ J2 {Qh^iS^nb^l 6^(aq,« + 

{D-SD)<\e^,\,\e^„\<D,a,P,a' 

X G';j-(e) &^(aq',a +aq.,5) G'J(e) ^^(aq'.fi + aq'.a)l(q>«)) 


E (QhHsh^ {ny\blc^,^blcl 

(D-5D)<\e^,\,\e^„\<D,a,P,a' 


E-{Ho + e, + e,,o) + “q'.^) 


(Oq'.a + Oqys)!^) 


E — {Hq -h £;, -f eqr,o) 




N 


q'./3 




E — 2eb 
1 


' ~ (”£ 9 ',! + £b + £g.o) 
t 1 

“qCoOqco p /, I , I , V 

^ ~ (,£g',0 + £b + £g,oj 


~ -2Q^S'^ 


p{D)a^\5D\ 1 


E-D-J E-2J 


Y^iQ^isy. 


(5) 


Adding the two contributions (taking A ~ 0), 


T(« 2Q'>S'>p(Dfy {5 + 2 (dT 7 ) } ^ 

( 6 ) 

Here we have taken E, Cq^ci ^ D, J and neglected them. 

If either the impurity is a S' = ^ spin, or the Kondo 
interaction is rotationally symmetric, the above contri¬ 
bution renormalizes the Kondo coupling constant. How¬ 
ever for S 4 with anisotropic coupling, new terms are 
generated and one needs to go to higher order diagrams 
to obtain the scaling of these new coupling terms. For 
S = 1/2 or for symmetric impurity coupling we thus have 

(7) 

Just as for the Kondo effect in graphene!^, owing to the 
change in the density of states with bandwidth (here 
p(e) = {\/2TTv'^p)\e\ = G|e|), we also need to consider the 
change in K due to the rescaling done in order to keep 
the total number of states fixed. This gives a contribu¬ 
tion K —^ K{D'/D), [D' = D — |(5iJ|). In addition, as 
we shall scale the bandwidth D to smaller values, the 
second term in Eq. Hmay be dropped. Thus 

SK ~ -2K^S{S + l)p{D)a^^ + 

= -K^ {2K'^a'^CDS{S + 1)/J -1). (8) 

Thus, as we decrease the bandwidth by integrating out 
the high energy excitations, the effective coupling K has 
an unstable fixed point at Kc = y^J/[2a^p(iJ)S'(S' -I- 1)]; 
or in other words, Kc ^ \JJ/S'^a^CD ^ J/S. Here we 
used D < J and G ^ l/(Ja)^. Clearly for K > Kc, 
the coupling flows to infinity independent of the nature 
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of coupling (ferromagnetic or antiferromagnetic), while 
for K < Kc, the coupling flows to zero. For anisotropic 
Kondo coupling we can show 


where a, j3 belong to the ground state manifold and 7 
belongs to excited state manifold. The eigenstates of the 
Kondo term are singlet |s) and triplet states |t, (0,±1)): 


5K. 


z,± 


-K. 


z,±- 


SD_ 

D 


2a^p{D)S{S 


Kl+K+K_ ' 
’ J 

(9) 


The two-parameter Kondo flow is therefore given by 


5K, 

5K± 


K± 


E± 

K± 


= const. 


( 10 ) 


A comparison of the Kondo effect in graphenei^, a 
bosonic spin batb^ and the Kitaev model are shown in 
Table H 


|s) = ^(|t,^)-U,^)), 

(17) 


(18) 

( ti-D-) + N-ilt)) : 

(19) 


( 20 ) 


Here f refers to the Kitaev spin and fi refers to the im¬ 
purity spin state. 


Antiferromagnetic Kondo coupling 


A. Stability of strong coupling point 


For the antiferromagnetic Kondo coupling case, ground 
state is the singlet state. As (s|F|s) = 0, i?i = 0 and 


The poor man’s scaling analysis is only valid for small 
Kondo couplings as the perturbation theory breaks down 
much before the critical value of the coupling. While we 
have shown that the coupling flows to larger values above 
the critical value Kc, it remains to be seen whether there 
is any other fixed point beyond Kc but less than the 00 . 
Below we study the model in the strong coupling limit 
and see if it is a stable fixed point. In the strong coupling 
limit, K is the largest energy scale and the impurity spins 
forms a singlet/triplet with the Kitaev spin at origin. 

We consider the Hamiltonian such that the Kondo 
term and the Kitaev model with one spin missing {Hk-) 
constitute the unperturbed Hamiltonian and Kitaev cou¬ 
pling to the site at origin is the perturbation: 


Ho = KS-ao + HK-, ( 11 ) 

l/ = J(aoX + aoM + aoV3^). (12) 


For antiferromagnetic Kondo coupling {K > 0), the 
ground state consists of a Kondo singlet of S and cto 
and the Kitaev model with one spin missing. The per¬ 
turbation term causes transitions from singlet to triplet 
states of the Kondo singlet. We use effective Hamiltonian 
schemers to include the effects of the perturbation terms 
within the projected ground state subspace. 

Hcs = e^^He-^^, (13) 


where Q is chosen such that the terms which take us out 
of the reduced Hilbert space are canceled order by order. 
This gives the reduced Hamiltonian as 


Heff = Ho + Hi+H 2 + 0{V^), 
{a\Hi\P) = {a\V\(3), 

(a|H2|/3) = i^(a|F|7)(7im 

7/“,/3 


1 

Ea — 


(14) 

(15) 


Ea — E^ 


(16) 


(s,K_|i?2|s,X^) = i ^(s,X-|F|t,A") 


t,K" 


X {t,Kl\V\s,K'_) 


1 


1 


Eq — Et E'o — Et 


■ ( 21 ) 


Here, AT- denotes the eigenstates of the Kitaev model 
with the spin at origin missing. Since change in energy 
of the Kitaev state is ^ J <C AT, we ignore their contri¬ 
bution in the energy denominators of the perturbation 
term. The matrix elements of H 2 are then 


{s,K_\H2\s,K'_) (22) 

^ ^(s|a“|t)(t|ao^|s)(if_|a“|i^")(if"|a^|A^) 

(23) 

- (1 “ \s){s\)cr^\s){K_\a2(7^\K'_) (24) 

o j2 

(25) 

Here, in cr“, the subscript a refers to a neighoring site 
of the origin in the direction of the a—bond. Thus, in 
the antiferromagnetic coupling case, the Kondo singlet 
decouples from the rest of the Kitaev model and a small 
interaction (J^/AT <C J) is generated between the Kitaev 
spin at the origin and the sins at the three neighbor¬ 
ing sites in the second order perturbation. The strong 
coupling fixed point is thus a stable fixed point and is 
equivalent to the Kitaev model with one site missing. 


Ferromagnetic Kondo coupling 

In the ferromagnetic Kondo coupling case, the triplet 
states form the ground state manifold. We perform de- 
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Graphene 

Z 2 bosonic spin bath with 
pseudogap density of states 
p(e) = C|e|. 

Kitaev, honeycomb lattice 

Kondo 

scaling 

Unstable intermediate 
coupling fixed pt. only for 
AFM coupling. Only AFM 
flows to strong coupling 
above unstable fixed pt. 

Flow direction is independent of 
the sign of magnetic impurity 
coupling. Unstable intermediate 
coupling fixed pt. for both FM 
and AFM. 

Scaling same as Z 2 bosonic 
spin bath case. However a 
topological transition is 
associated with the 
unstable fixed point. 


Table I. Comparison of Kondo effect in graphene, a Zi bosonic spin bath with a pseudogap density of states and the Kitaev 
model on the honeycomb lattice. 


generate perturbation theory to get the effective Hamil¬ 
tonian: 


= Jl4! ,K'_\V%K.) (26) 

= Y.^t'\oZ\t){K'_\al\K.). (27) 


If we calculate the matrix elements of (f'|crQ |t), these ma¬ 
trices are just the spin —1 matrices: 


/ 0 ^ 0 



-rV — 


/ 0 o_ 


V2 


0 


0 0 


1 0 0 

CTg = ( 0 0 0 

0 0-1 


and the Hamiltonian in the reduced subspace becomes 
i7i = (28) 

where Sg represents the spin -1 at the origin. 

Thus for ferromagnetic impurity coupling, the new 
terms which couple the triplet and the rest of the Kitaev 
model are similar to the original Kitaev coupling and of 
the same strength. We get a Kitaev-like model with a 
spin—1 at the origin and spin—1/2 elsewhere. Here the 
Kondo triplet does not decouple from rest of the Kitaev 
model in the strong coupling limit and does not lend itself 
to a simple treatment, unlike the corresponding antifer¬ 
romagnetic case. 


B. Topological transition 

A remarkable property of the Kondo effect in Kitaev 
model is that the unstable fixed point is associated with 
a topological transition from the zero flux state to a finite 
flux state. The strong antiferromagnetic coupling limit 


amounts to studying the Kitaev model with a missing site 
or cutting the three bonds linking this site to the neigh¬ 
bors. It was shown in Kitaev’s original paper— that such 
states with an odd number of cuts are associated with 
a finite flux, and also that these vortices are associated 
with unpaired Majorana fermions and have non-abelian 
statistics under exchange. It has also been shown nu¬ 
merically for the gapless phasei^ that the ground state 
of the Kitaev model with one spin missing has a finite 
flux pinned to the defect site. We argue the existence 
of a localized zero energy Majorana mode from the de¬ 
generacy of the ground state in presence of impurity spin 
and elucidate on the nature of this zero mode. 

For the Hamiltonian H — Hq + Vk, the three plaque- 
ttes Wi, W 2 and W 3 (Fig. [3]) that touch the impurity 
site are no longer associated with conserved flux oper¬ 
ators, while the flux operators that do not include the 
origin remain conserved. The three plaquette operator 
Wo = W 1 W 2 W 3 is still conserved and Wq = 1 in the 
ground state of the unperturbed Kitaev model. 

We now define composite operators = W 2 W 3 S''^, 
ry = WsWiSy and = Wi WsS”^ ( 5“ are the Pauli spin 
matrices corresponding to the impurity). Remarkably, 
these composite operators represent conserved quantities 
for arbitrary values of the impurity coupling. The r“’s 
do not commute with each other and instead obey an 
SU{2) algebra, [r“,r^] = 2i€ai3-yT'^ ■ This SU{2) sym¬ 
metry, which is exact for all couplings is realized in the 

spin-1/2 representation ^(t“)^ = 1^ Clearly, all eigen¬ 
states, including the ground state are doubly degenerate 
(corresponding to = ± 1 ), and this applies also to the 
strong coupling limit. 

In the strong antiferromagnetic coupling limit Jk 
00 , the low energy states will be the ones in which the 
spin at the origin forms a singlet |0) with the impurity 
spin, !'(/)) = \%pK~) ® |0). Here \'il)K~) represents the low 
energy states of the Kitaev model with the spin at the 
origin removed. To see the action of the SU{2) symme¬ 
try generators on these states, we note that they can be 
written as r“ = W“ C) erg 0 S'“ and W“ do not involve the 
components of the spin at the origin, dg . We then have 
t“|' 0) = —(W“|'07f“)) 0 | 0 ). So, in the strong coupling 
limit, the symmetry generators act non-trivially only in 
the Kitaev model sector, implying that the low energy 
states of the Kitaev model with one spin removed are all 
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Figure 3. Schematic of the three unpaired b-Majorana 
fermions formed as a result of cutting the links to the Ki- 
taev spin at the origin. Any two of the three can be given an 
expectation value (dotted bond). 


doubly degenerate, with the double degeneracy emerging 
from the Kitaev sector. This is also true for the zero- 
energy mode in the single particle spectrum: the two 
degenerate states correspond to the zero mode being oc¬ 
cupied or unoccupied. 

Let us examine the structure of the zero mode. Re¬ 
moving a Kitaev spin creates three unpaired 6 —Majorana 
fermions at the neighboring sites, say, 63 , bf and (Fig. 
[3|). Now ibfb 2 is conserved and commutes with all the 
conserved flux operators Wi but not with the two other 
combinations ib^b^ and ib^bf. So, we can choose a gauge 
where the expectation value of ibfb^ is equal to -|-1 such 
that these two 6 —modes drop out of the problem and 
we equivalently have one unpaired 6 —Majorana fermion. 
The unpaired Majorana has dimension and there¬ 
fore, there must be an unpaired Majorana mode in the 
c—sector (again of dimension v^) so that together these 
two give the full (doubly degenerate) zero energy mode. 
Also, while the 63 mode is sharply localized, the wave 
function of the c mode can be spread out in the lattice. 

For the ferromagnetic case, while the strong coupling 
limit also leads to a model with doubly-degenerate levels, 
we are unable to explicitly identify a zero energy unpaired 
Majorana fermion and not address the question as to 
whether a nontrivial Z 2 flux can be associated with closed 
paths enclosing the S' = 1 defect. For this purpoose, 
we perform a numerical exact diagonalization analysis 
below. 


IV. NUMERICAL STUDIES 

We have used a modified Lanczos algorithm to calcu¬ 
late the ground state properties of a finite Kitaev frag¬ 
ment exchange-coupled to an external S = 1/2 impurity 
spin as discussed in the previous sections. 


For the antiferromagnetic case, where we already know 
that the strong-coupling limit K/J —>• 00 is associated 
with a nontrivial flux Wq = —1 at the defect site, 
an exact-diagonalization calculation with a fragment as 
small as three hexagons (open boundary conditions, im¬ 
purity spin coupled to central site) is sufficient to con¬ 
firm Wq = — 1. Figure m shows the expectation value of 
the flux operator W and total spin Srot = S + Si as & 
function of the impurity coupling K. It is seen that the 
Wo = 1 for AT = 0 (i.e in the pure Kitaev case) as it 
should be, since the ground state of Kitaev model is flux 
free, whereas it changes to -1 as K is increased, implying 
a finite Z 2 flux at the origin. Within numerical accu¬ 
racy, it also appears that the topological transition from 
Wo = 1 to Wo = — 1 practically occurs at the same value 
of K at which a bound singlet state is formed between 
the impurity spin and the Kitaev host. For negative val¬ 
ues of K for this three-hexagon fragment. Wo appears to 
stay close to one implying a flux-free state even as the 
total spin at the defect site begins approaching S = 1. 
We suspect this anomalous result is not generally true 
for larger fragments and may have originated from the 
presence of a large number of boundary spins that in¬ 
teract with the central spin. We thus performed exact- 
diagonalization calculations with a larger fragment with 
six hexagons (open boundary conditions, impurity spin 
coupled to central site). As we increase the ferromag¬ 
netic impurity coupling, we clearly observe the Wq = 1 
to Wo = — 1 topological transition. Once again, the topo¬ 
logical phase transition and the magnetic transition (at 
which the total spin at the defect site becomes S' = 1, 
are practically coincident. 


V. TWO-CHANNEL KONDO BEHAVIOR 

In Ref.[ll the temperature and magnetic field depen¬ 
dences of the magnetic susceptibility of the Kitaev model 
with a missing site were obtained as Ximp(F) ^ ln(l/T) 
and Ximp(/i) lii(l//i), respectively, which bears strik¬ 
ing resemblance to the low-temperature impurity suscep¬ 
tibility in the two-channel Kondo model^i^. This is 
not a mere coincidence. In the two-channel Kondo prob¬ 
lem, it is long knowni^ that the low energy physics is 
described by a model of a localized, zero energy Majo¬ 
rana fermion interacting with a band of dispersing Majo¬ 
rana fermions with a finite density of states at the Fermi 
energy. Likewise, in the Kitaev model with a missing 
site, we can clearly identify a localized zero energy Ma¬ 
jorana 6 —fermion coexisting with dispersing Majorana 
c—fermions with a nonvanishing density of states^^ at 
zero energy. A nonvanishing density of states for the 
dispersing Majoranas is a very unusual result for a hon¬ 
eycomb lattice, and is associated with the fact that the 
missing site is associated with a finite Z 2 flux. If one 
instead estimates the density of states for the dispers¬ 
ing fermions in the absence of a tt— flux, the density of 
states would vanisbi^ at the Fermi energy; indeed, this 
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Figure 4. (Color online) (Top) Evolution of expectation value 
of the flux operator (Wo) and the total spin at the impu¬ 
rity site with K evaluated on a three hexagon fragment. The 
impurity spin is coupled to the central spin and the flux is 
calculated over the boundary. The inset shows the slight de¬ 
crease of Wo from 1 for ferromagnetic coupling whereas for 
antiferromagnetic coupling Wo changes to -1. (Bottom) Wo 
and total spin at the impurity site evaluated for a six hexagon 
fragment showing the transition to Wo = -1 state for ferro¬ 
magnetic coupling. 


would be the case in graphene. To compute the magnetic 
susceptibility, we choose a gauge where is the zero en¬ 
ergy localized Majorana fermion. Using = ic^b^ the 
magnetic susceptibility of the defect can be expressed as 

Ximp = T^^-^^ln(l/r), (29) 

^ IVn - €k ll'n 

n,k 


ln(l/h). For the ground state entropy, one notes that the 
dispersing Majorana fermions have zero entropy at T = 0 
while the localized zero energy Majorana fermion has a 
finite entropy S = (1/2) ln(2). This again agrees with the 
two-channel Kondo result 

VI. DISCUSSION 


In summary, we showed that the problem of spinless 
and spinful defects in the honeycomb S' = 1/2 Kitaev 
model can be approached from a more general “Kondo 
perspective” of local exchange coupling of external para¬ 
magnetic impurities with a host spin. On one hand, such 
an approach gives us a new class of Kondo effects where 
the magnetic binding-unbinding transition is accompa¬ 
nied by a change of topology of the ground state. On 
the other hand, some intriguing recent observations, such 
as logarithmic singularities in the magnetic response of 
Kitaev models with vacancies, are now recognizable as 
familiar Kondo stories - in this case, we note a remark¬ 
able similarity with the two-channel Kondo problem. It 
would be interesting to study a lattice of vacancies in 
the Kitaev model from the perspective of a two-channel 
Kondo lattice. One would like to better understand the 
Kitaev model with a S' = 1 defect. This nonintegrable na¬ 
ture of this problem prevents us from repeating the kind 
of analysis one could make for the vacancy case where 
similarity with the two-channel Kondo problem was ob¬ 
served. Our numerical approach, based on exact diago- 
nalization calculations of relatively small fragments, can¬ 
not answer questions such as the density of states of low 
energy excitations. Another direction for future study 
would be to consider a more general Kitaev-Heisenberg 
model and track the Kondo effect as a function of the rel¬ 
ative strengths of Kitaev and Kondo interactions. This 
should give insights into magnetic impurity response in 
Kitaev candidate materials where Kitaev and Heisenberg 
interactions are believed to compete with each other. 
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